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We present a measurement of the fractions Fo and F+ of longitudinally polarized and right-handed 
W bosons in top quark decays using data collected with the CDF II detector. The data set used 
in the analysis corresponds to an integrated luminosity of approximately 318 pb -1 . We select ti 
candidate events with one lepton, at least four jets, and missing transverse energy. Our helicity 
measurement uses the decay angle 8*, which is defined as the angle between the momentum of the 
charged lepton in the W boson rest frame and the W momentum in the top quark rest frame. The 
cos 9* distribution in the data is determined by full kinematic reconstruction of the ti candidates. 
We find F = 0.85±°;22(stat) ±0.06 (syst) and F+ = 0.05±° \H (stat) ±0.03 (syst), which is consistent 
with the standard model prediction. We set an upper limit on the fraction of right-handed W bosons 
of F + < 0.26 at the 95% confidence level. 
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I. INTRODUCTION 

In 1995 the top quark was discovered at the Tevatron 
proton-antiproton collider at Fermilab by the CDF and 
D0 collaborations [1, 2]. It is the most massive known el- 
ementary particle and its mass is currently measured with 
a precision of about 1.3% [3, 4]. However, the measure- 
ments of other top quark properties are still statistically 
limited, so the question remains whether the standard 
model successfully predicts these properties. This paper 
addresses one interesting aspect of top quark decay, the 
hclicity of the W boson produced in the decay t — ► W+ b. 

At the Tevatron collider, with a center-of-mass energy 
t/s = 1.96 TeV, most top quarks are pair-produced via 
the strong interaction. In the standard model the top 
quark decays predominantly into a W boson and a b 
quark, with a branching ratio close to 100%. The V — A 
structure of the weak interaction of the standard model 
predicts that the W+ bosons from the top quark decay 
t — ► W+ b are dominantly either longitudinally polarized 
or left-handed, while right-handed W bosons are heavily 
suppressed and are forbidden in the limit of massless b 
quarks. 

As a consequence of the Goldstone boson equivalence 
theorem [5, 6], the decay amplitude to longitudinal W 
bosons is proportional to the Yukawa coupling of the top 
quark; therefore, the decay rate scales with mf , where m t 
is the top quark mass. The longitudinal decay mode of 
the W boson is thereby linked to the spontaneous break- 
ing of the electroweak gauge symmetry. The decay rate 
to transverse W bosons is governed by the gauge coupling 
and increases only linearly with m t [7]. The fraction of 
longitudinally polarized W bosons is defined by 



F 



T(t W+b) 



T(t -» W+b) + T(t -» W+b) + r(t -» W+b) 



(1) 



where W+ stands for a longitudinally polarized W+ bo- 
son, W+ for a left-handed W+ boson, and W+^ for a 
right-handed W+ boson. The corresponding definitions 
for the W~ boson are implied. In leading-order pcrtur- 

2 

bation theory F a is predicted to be F a = " lt . 2 [8], 

where mw is the mass of the W boson. Using mw — 
80.43 GeV/c 2 [9] andm t = (172.5±2.3) GeV/c 2 [3], gives 
F = 0.697 ± 0.007, where the given uncertainty is only 
due to the uncertainty in the top quark mass. Next-to- 
lcading-order corrections to the total decay width and the 
partial decay width into longitudinal W bosons amount 
to about -10% [10-19]. However, the fraction of longitu- 
dinal W bosons is only negligibly changed. 
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A significant deviation from the predicted value for Fq 
or a nonzero value for the right-handed fraction F + could 
indicate new physics. Left-right symmetric models [20], 
for example, lead to a significant right-handed fraction 
of W bosons in top quark decays. Such a right-handed 
component (V + A coupling) would lead to a smaller left- 
handed fraction, while the longitudinal fraction Fo would 
change insignificantly. Since the decay rate to longitudi- 
nal W bosons depends on the Yukawa coupling of the 
top quarks, the measurement of F is sensitive to the 
mechanism of electroweak symmetry breaking. Alterna- 
tive models for electroweak symmetry breaking, such as 
topcolor-assisted technicolor models, can lead to an al- 
tered F fraction [21, 22]. 

The W boson polarization manifests itself in the decay 
W — > Iv in the angle 8*, which is defined as the angle 
between the momentum of the charged lepton in the W 
rest frame and the W momentum in the top quark rest 
frame. For a longitudinal fraction F , a right-handed 
fraction F + , and a left-handed fraction F_ = 1 — F + — F , 
the cos 9* distribution is given by [8] : 



<IX - (i-f+-fb)-|(i-cosr) 2 



dcosd* 



+ (F )- i (l-cos 2 r) 
+ (F + )-|(l + cosr) 2 . 



(2) 



In this analysis, the W hclicity fractions are measured 
in a selected sample rich in tt events where one lepton, 
at least four jets, and missing transverse energy are re- 
quired [23]. In order to calculate 9*, all kinematic quan- 
tities describing the tt decays have to be determined. 

Previous CDF measurements of the W helicity frac- 
tions in top quark decays used either the square of the 
invariant mass of the charged lepton and the b quark 
jet [24-26] or the lepton pt distribution [25, 27] as 
a discriminant. The D0 collaboration used a matrix- 
element method to extract a value for F [28]; in a sec- 
ond analysis the reconstructed distribution of cos 9* [29] 
was utilized to measure F+. The previous measurement 
by CDF was F = 0.74+_oil [25], while D0 measured 
F = 0.56 ± 0.31 [28]. The CDF collaboration also mea- 
sured the current best upper limit of F + < 0.09 at the 
95% confidence level [26]. 

The organization of this paper is as follows. Section 
II describes the detector system relevant to this analysis. 
Section III illustrates the event selection of the tt candi- 
dates. The signal simulation and background estimation 
are given in Section IV. In Section V we describe our 
method to fully reconstruct tt pairs. The extraction of 
the helicity fractions is presented in Section VI. Section 
VII discusses the systematic uncertainties. Finally, the 
results and conclusions are given in Section VIII. 
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II. THE CDF II DETECTOR 

A detailed description of the Collider Detector at Fer- 
milab (CDF) can be found elsewhere [30]. A coordinate 
system with the z axis along the proton beam, azimuthal 
angle <j>, and polar angle 9 is used. The azimuthal angle 
is defined with respect to the outgoing radial direction 
and the polar angle is defined with respect to the pro- 
ton beam direction. The transverse energy of a particle 
is defined as Et — EsinO. Throughout this paper we 
use pseudorapidity defined as r\ = — ln(tan|). The pri- 
mary detector components relevant to this analysis are 
those which measure the energies and directions of jets, 
electrons, and muons and are briefly described below. 

An open-cell drift chamber, the central outer tracker 
(COT) [31], and a silicon tracking system arc used to 
measure the momenta of charged particles. The CDF II 
silicon tracker consists of three subdetectors: a layer of 
single-sided silicon microstrip detectors [32] glued on the 
beam pipe, a five layer double-sided silicon microstrip de- 
tector (SVX II) [33], and intermediate silicon layers [34] 
located at radii between 19 and 29 cm which provide link- 
ing between track segments in the COT and the SVX II. 
In the analysis presented in this article, the silicon tracker 
is used to identify jets originating from b quarks by re- 
constructing secondary vertices. The tracking detectors 
are located within a 1.4 T solenoid. Electromagnetic and 
hadronic sampling calorimeters [35-37], which have an 
angular coverage of |7y| < 3.6, surround the tracking sys- 
tem and measure the energy flow of interacting particles. 
They are segmented into projective towers, each one cov- 
ering a small range in pseudorapidity and azimuth. For 
electron identification the electromagnetic calorimeters 
are used, while jets are identified through the energy they 
deposit in the electromagnetic and hadronic calorimeter 
towers. The muon system [38] is located outside of the 
calorimeters and provides muon detection in the range 
| t? | < 1.5. Muons penetrating the five absorption lengths 
of the calorimeters are detected in planes of multi-wire 
drift chambers. Since the collision rate exceeds the tape 
writing speed by five orders of magnitude, CDF has a 
three-level trigger system which reduces the event rate 
from 1.7 MHz to 60 Hz for permanent storage. The first 
two levels of trigger are implemented by special-purpose 
hardware, whereas the third one is implemented by soft- 
ware running on a computer farm. 

III. SELECTION OF ti CANDIDATE EVENTS 

In the decay channel considered in this analysis, one 
top quark decays semileptonically and the second top 
quark decays hadronically, leading to a signature of one 
charged lepton, missing transverse energy resulting from 
the undetected neutrino, and at least four jets. Can- 
didate events are selected with high-p T lepton triggers. 
The electron trigger requires a COT track matched to an 
energy cluster in the central electromagnetic calorime- 



ter with Et > 18 GeV. The muon trigger requires a 
COT track with pt > 18 GeV/c matched to a track 
segment in the muon chambers. After offline recon- 
struction, we require exactly one isolated electron can- 
didate with Et > 20GcV and \r/\ < 1.1 or exactly 
one isolated muon candidate with pt > 20GeV/c and 
|r/| < 1.0. An electron or muon candidate is considered 
isolated i f the Et not as signed to the lepton in a cone 
of R = \J (Ar)) 2 + (A(f>) 2 =0.4 centered around the lep- 
ton is less than 10% of the lepton E T or p T , respectively. 
Jets are reconstructed by summing calorimeter energy 
in a cone of radius R = 0.4. The energy of the jets 
is corrected [39] for the r\ dependence of the calorime- 
ter response, the time dependence of the calorimeter re- 
sponse, and the extra deposition of energy due to mul- 
tiple interactions. Candidate jets must have corrected 
Et > 15 GeV and detector \r]\ < 2.0. Detector rj is 
defined as the pseudorapidity of the jet calculated with 
respect to the center of the detector. Events with at least 
four jets are accepted. At least one of the jets must be 
tagged as a 6-jet by requiring a displaced secondary ver- 
tex within the jet [40]. The missing Et (^?t) is defined 

by 

$T = -J2 E T^ (3) 

i 

i = calorimeter tower number with |?7| < 3.6, 

where hi is a unit vector perpendicular to the beam axis 
and pointing at the i th calorimeter tower. We also de- 
fine E T = | $t\- Because this calculation is based on 
calorimeter towers, Et has to be adjusted for the effect 
of the jet corrections for all jets with Et > 8 GeV and 
detector |ry| < 2.5. In events with muons, the transverse 
momentum of the muon is added to the sum, and a cor- 
rection is applied to remove the average ionization energy 
released by the muon in traversing the calorimeter. We 
require the corrected Et to be greater than 20 GeV. 

Additional requirements reduce the contamination 
from background. Electron events are rejected if the 
electron stems from a conversion of a photon. Cosmic 
ray muon events are also excluded. To remove Z boson 
events, we reject events in which the charged lepton can 
be paired with any more loosely defined jet or lepton to 
form an invariant mass consistent with the Z peak, de- 
fined as the range 76 GeV/c 2 to 106 GeV/c 2 . After these 
selection requirements we find 82 ti candidates in the se- 
lected sample corresponding to an integrated luminosity 
of 318 pb" 1 . 



IV. SIGNAL SIMULATION AND 
BACKGROUND ESTIMATION 

In order to determine the resolution of the kinematic 
quantities of the reconstructed ti pair, as well as to de- 
termine certain background rates, we utilize Monte Carlo 
simulations. The generated events are passed through 
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the CDF detector simulation [41] and are reconstructed 
in the same way as the measured data. 

The simulated ti signal sample was generated with the 
pythia generator [42] using a top quark mass of m t — 
178 GeV/c 2 which was the world average [43] of Run I. 
The values of Fq and F + used in our standard model 
simulation are 0.7 and 0.0 respectively. To check the 
assumption that neither the efficiency nor the resolution 
due to the reconstruction depend on the values of Fq and 
F+, we use a customized version of the herwig Monte 
Carlo program [44] in which the hclicity of one W boson 
is fixed to be longitudinal, left-handed, or right-handed. 

The selected ti candidate sample contains a certain 
level of background contamination. Among the 82 ob- 
served events, we predict a background of 10.3 ±1.9 
events [23]. The dominant sources are W production 
in association with a quark-antiquark pair (31%), e.g. 
qq' — > Wgg with g — > bb (cc) and g — > q"q", "mistagged" 
events (24%), in which a jet is erroneously tagged as 
a fe-jct, and events where no W boson (non-W events) 
is produced (36%), e.g. direct bb production with addi- 
tional gluon radiation. Additional sources are diboson 
(WW, WZ, ZZ) production (4.5%) and single-top pro- 
duction (4.5%). The non-FT and mistag fractions arc 
estimated using lepton trigger data. The W plus heavy 
flavor fraction is extracted using a sample of events simu- 
lated with ALPGEN [45]. The diboson and single-top rates 
are predicted based on their theoretical cross sections [46] 
and acceptances and efficiencies, which are derived from 
PYTHIA and MADEVENT [47] simulations. 



V. FULL RECONSTRUCTION OF tt PAIRS 



the semileptonically decaying top quark, we consider all 
combinations of the four-momentum of one of the se- 
lected jets and the four-momentum of the W boson. The 
hadronically decaying W boson is then reconstructed by 
combining the four-momenta of two of the selected jets 
not assigned to the semileptonically decaying top quark. 
Adding the four-momenta of this W boson and of one 
of the remaining jets results in the hadronically decaying 
top quark. This procedure leads to \ • Nj ets \ / (Nj ets — 4)! 
different hypotheses for each event. 

For simulated events it is possible to determine the 
hypothesis which is closest to the true event. This "best 
hypothesis" is defined as the hypothesis for which the 
deviation of the reconstructed top quarks and W bosons 
from the generated particles in the rj-cj) plane is minimal. 
Since this is not possible for measured data, we determine 
for each hypothesis a quantity if which gives a quanti- 
tative estimate of how well this hypothesis matches the 
ti pair assumption, and we choose the hypothesis with 
the highest value of \&. Constraints on the mass of the 
hadronically decaying W boson, on the mass difference 
between both reconstructed top quark masses, on how 
Mike the jets assigned as 6-jets in the t — ► Wb decays 
are, and on the reconstructed Et of the two top quarks 
enter the computation of \F 

We define W as 



•ft , 



(4) 



where Je is the sum of the transverse energies of the 
two top quarks divided by the total Et of the event 
including Jfi^: 



The measurement of cos 9* is based on fully recon- 
structing the top quarks through the four-momenta of the 
decay products. The challenge for the full reconstruction 
is to assign the observed jets to the decay products of the 
hadronically decaying W boson or the jets resulting from 
the b quarks from the top-quark decays. All possible as- 
signments have to be considered. Thus, in each event 
there exist numerous hypotheses for the reconstruction 
of the tt pair. At the top quark reconstruction level, ex- 
tra jet corrections are applied. The calorimeter energy is 
corrected to correspond to the energy of the traversing 
particle, the underlying event energy is subtracted, and, 
finally, the energy that is radiated outside the jet cone is 
added. The pt vector of the neutrino is derived from • 
To calculate the z-component of the neutrino momentum, 
a quadratic constraint using the W ' — > iv decay kinemat- 
ics is used, with the assumption that the W boson mass 
equals the pole mass of 80.43 GeV/c 2 . If the solution 
of the equation is complex, the real part of the solution 
is taken; otherwise the solution with the smaller value 
of \p z ,v\ is used. Adding the resulting four- momentum 
of the neutrino and the four-momentum of the charged 
lepton leads to the correct W boson four-vector in 78% 
of simulated events. In order to get all hypotheses for 



^PT,jet+ + Et,£ 

where PT,t->6fc and PT.t^bjj are the reconstructed trans- 
verse momenta of the semileptonically and hadronically 
decaying top quarks and mt^beu and m t ^bjj are the re- 
spective reconstructed top quark masses. The quantity 
SpT.jet is the sum of the transverse momenta of the four 
jets used in the ti event hypothesis. The transverse en- 
ergy of the charged lepton is indicated with i?T,<- The 
motivation for the definition of Je is that the Et of the 
top quarks is approximately equal to the Et of the en- 
tire event. The mean value Je of the Je distribution, 
obtained from the best hypothesis for each event of a ti 
Monte Carlo simulation, is determined to be 1.014. 
The quantity y 2 is defined as 



X 



[m W - 



2 -i- 2 ' w 



Ami 



where mw^jj is the reconstructed mass of the hadron- 
ically decaying W boson and m t ^biv and m t ^bjj are 
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TABLE I: Percentage of tt events that are reconstructed 
within a particular ^2 AR, as defined in Eq. 8. 



fraction [%] 


best 


30.2 


X>i?< 1.5 


41.5 


J2 A R< 3-0 


57.9 


^ AT? < 4.5 


66.4 



the reconstructed mass of the semileptonically decaying 
top quark and the hadronically decaying top quark, re- 
spectively. The reconstructed mass of the hadronically 
decaying W boson should be equal to the mean value 
m w ^jj within the resolution o. mw ^ ij and the difference 
between both top quark masses should be zero within 
the resolution a/\„ lt . The values fhw^jj = 79.5 GeV/c 2 , 
a mw ^ j3 = 10.2GcV/c 2 , and a Amt = 30.3GeV/c 2 that 
we use are obtained from the corresponding mass dis- 
tributions using the best hypothesis of fully simulated 
tt events. The mass resolutions are dominated by the 
uncertainties in the jet energy reconstruction. The jet 
energy scale is determined from dijet data events and 
simulated samples and checked using 7+jet and Z+jet 
events [4, 48]. The value for rhw^jj deviates from the 
measured W boson pole mass mw — 80.43 GeV/c 2 . The 
deviation is within the systematic uncertainties of the 
applied jet corrections. 

The quantity is a measure of how 6-like the two jets 
assigned as such by the event reconstruction are, and is 
defined as: 

P b = (-logPt^-logPt^O-lO^ , (7) 

where Vt^biv and Vt^bjj are the probabilities that the 
jets chosen to be the 6-jets from the semileptonically and 
hadronically decaying top quark are consistent with the 
hypothesis of a light quark jet with zero lifetime. This 
probability is calculated from the impact parameter of 
the tracks assigned to the jet in the r-<fi plane [49]. The 
negative logarithm of that probability leads to large val- 
ues for 6-jets and small values for light flavor jets. How- 
ever, since a reconstructed secondary vertex is a stronger 
indication for 6-jets than the probability based on the 
impact parameter, the quantity Pb should be given a 
higher weight when there are secondary vertex tagged 
jets. Since — \ogV nearly always takes values smaller 
than 10, the logarithmic sum is multiplied by the factor 
10 JVta «, where N tag is the number of 6-tagged jets (either 
0, 1 or 2). 

In order to estimate the quality of the criterion 
for choosing the most probable event reconstruction 
based on the quantity Monte Carlo studies are per- 
formed. We examine the sum of the distances in the 
rj-tfi plane associated with the semileptonically decay- 
ing top quark {AR t ^uv), the hadronically decaying top 
quark (AR t ^ b jj), and the hadronically decaying W bo- 
son (ARw^jj)- 



• data 




FIG. 1: Measured cos#* ec distribution shown together with 
the estimated signal and background. The tt signal events 
are modeled by the standard model Monte Carlo generator 
pythia with m t = 178 GeV/c 2 . 



^AR = AR t ^ Miy + AR t ^ bj3 + ARw^jj. (8) 

The distance AR between a generated ("gen") and 
a reconstructed ("rec") particle is given by AR = 
V ('/'gen - 4>rcc) 2 + (%cn — ?7rcc) 2 - Table I shows how of- 
ten our selected hypothesis has a value of AR below a 
given value. We also state the fraction of events in which 
the chosen hypothesis is the "best hypothesis" which is 
defined for each event as the hypothesis with the smallest 
value of A-R- 

Our reconstruction method yields cos 6** resolutions 
comparable to other methods used in previous CDF mea- 
surements [48]. In addition the present approach allows 
the inclusion of events with more than four jets in a con- 
sistent way. 

Figure 1 shows the distribution of the measured cos 9* 
compared to the estimated signal and background distri- 
butions. 



VI. EXTRACTION OF F AND F+ AND 
DETERMINATION OF THE DIFFERENTIAL it 
PRODUCTION CROSS SECTION 

Since the number of events in the data set is small, we 
do not simultaneously extract the fraction of longitudi- 
nally polarized and right-handed W bosons. We either 
fix F + to and fit for Fq, or we fix Fq to its expected 
value and fit for F + . Thus, only one free parameter is 
used in each fit. 

To extract the single free parameter [Fq or F+), we 
use a binned maximum likelihood method. The expected 
number of events in each bin is the sum of the expected 
background and signal. The latter is calculated from 
the theoretical cos 8* distributions (Eq. 2) for the three 
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hclicitics of the W boson. Integrating Eq. 2 for each bin 
i separately leads to a linear dependence of the expected 
number of signal events /x s on F and F + : 

tf* cx (1 - F - F + ) ■ fr + (F ) ■ /° + (F + ) ■ /+. (9) 
Here /°, /~ and /+ are defined as: 

rb, o 

fi = -{1 - cos 2 9*)d cos 9*, (10) 

fi = j^\{l-cos9*) 2 dcos9\ (11) 

ft = ^ hi + cosO*) 2 d cos 9*, (12) 

Jo, 8 

where <Xj (6j) is the lower (upper) edge of the i th bin. 

As mentioned above, the reconstruction of the ti pro- 
cess is not perfectly efficient. Thus, in order to calculate 
the number of signal events ^ si e> obs expected to be ob- 
served in a certain bin after the reconstruction, we con- 
sider acceptance and migration effects: 

^^^■ei-S^k). (13) 

i 

The migration matrix element S(i, k) gives the proba- 
bility for an event which was generated in bin i to occur 
in bin k of the reconstructed cos#* distribution. Since 
the acceptance depends on cos#*, we weight the contri- 
bution of each bin i with the efficiency ej. Both e, and 
S(i, k) are determined using the standard model Monte 
Carlo generator pythia, assuming that ej and S(i, k) are 
independent of Fq and F + . This assumption has been 
verified using the customized herwig samples described 
above, which have fixed W helicities. 

With the number of expected events and the number 
of observed events in each bin, we minimize the negative 
logarithm of the likelihood function by varying the free 
parameter F or F + . 

In addition, an upper limit for F + at the 95% confi- 
dence level (CL) is computed by integrating the likeli- 
hood function L(F + ). Since a Bayesian approach is pur- 
sued, we integrate only in the physical region < F + < 
0.3 applying a prior distribution which is 1 in the interval 
[0,0.3] and elsewhere. 

In order to compare our observations with theory, the 
background estimate is subtracted from the selected sam- 
ple. To correct for acceptance and reconstruction effects, 
a transfer function r is calculated. The value r, for bin i 
is the ratio of the normalized number of theoretically ex- 
pected events and the normalized number of events after 
applying all selection cuts and performing the reconstruc- 
tion. For this calculation we use the fit result of Fo or 
F+. Multiplying the background-subtracted number of 
events in bin i with leads to the unfolded distribu- 
tion. Subsequently, this distribution is normalized to the 
ti production cross section of cr t j = 6.1 ± 0.9 pb [50, 51] 
assuming m t = 178 GeV/c 2 , which yields the desired 
distribution of the differential cross section. 



VII. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties caused by theoretical 
modeling, detector effects, and the analysis method have 
been studied using ensembles of simulated data samples. 
Each sample is made up of signal and background events 
drawn from the respective templates. The values for Fq 
and F + are extracted using the same method as for the 
observed data sample. The systematic uncertainty for a 
certain source is then given by comparing the mean of 
the resulting Fo and F + distributions of the correspond- 
ing ensemble with the default values. 

We account for possible bias from Monte Carlo model- 
ing of ti events by comparing herwig and pythia event 
generators. 

The contribution of the parton distribution function 
(PDF) uncertainty is determined by re-weighting the ti 
events generated with CTEQ5L [52] for different sets of 
PDFs. We add in quadrature the difference between 
MRST72 and MRST75 [53] and between the 20 pairs of 
CTEQ6M eigenvectors. 

To estimate the influence of initial-state and final-state 
radiation, we use templates from pythia Monte Carlo 
simulations in which the parameters for gluon radiation 
are varied to produce either less or more initial or final- 
state radiation [48] compared to the standard setup. The 
uncertainty due to the jet energy scale is quantified by 
varying the jet energy scale within its ±lcr uncertain- 
ties [39]. We also investigate whether our method to 
choose one hypothesis for each single event contributes 
significantly to the total uncertainty. Since the probable 
influence due to the x 2 and /b terms in the computation 
of the quantity ^ is already considered by varying the jet 
energy scale, we study the impact of Pb by omitting this 
term. To estimate the contribution of the background 
rate uncertainty, we simultaneously add or subtract, re- 
spectively, the values of one standard deviation of the es- 
timated rates for the different processes. The uncertainty 
due to the background shape uncertainty is estimated by 
using each shape of the dominant three background dis- 
tributions alone instead of using a composite of these 
shapes. 

The uncertainties are listed in Table II. The largest 
contribution to the systematic uncertainty arises from the 
jet energy-scale uncertainty, followed by the uncertainty 
on the background shape. 

Since the fraction of longitudinally polarized W bosons 
depends explicitly on the top quark mass, we do not 
include this dependence into the systematic uncertain- 
ties, but present our measurement assuming a certain 
top mass, namely 178 GeV/c 2 . 

However, we investigate the dependence of the mea- 
sured F and F + on the top quark mass. For a shift of 
+5 GeV/c 2 (—5 GeV/c 2 ) in the top quark mass we esti- 
mate a deviation in F of +0.017±0.007 (-0.017±0.007), 
which corresponds within the errors to the theoretical 

2 

prediction F = 2m "^ m 2 ■ The standard model predicts 
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FIG. 2: Extraction of the longitudinal (Jo) and right-handed (F+) fraction. For both fits Fo and F+ are used as single free 
parameter. In each case the other parameter is set to its expected standard model value. a,b) Negative log likelihood as 
a function of Fo or F+. c,d) Binned cos 9* distribution for data, corrected for acceptance and reconstruction effects. The 
distributions corresponding to the fit results Fo = 0.85 and F + = 0.05 are shown as a continuous functions. The dashed curve 
shows the theoretical prediction for Fo = 0.7 or for F + = 0.3. 



TABLE II: Summary of systematic uncertainties. The total 
uncertainty is calculated by adding all the individual uncer- 
tainties in quadrature. 



Uncertainties 



Source 


- AFo 


+ AF„ 


- AF+ 


+ AF+ 


Monte Carlo gen. 


0.022 


0.022 


0.010 


0.010 


Parton distribution functions 


0.017 


0.017 


0.006 


0.006 


Initial-state radiation 


0.010 


0.010 


0.007 


0.007 


Final-state radiation 


0.005 


0.005 


0.002 


0.002 


Jet energy scale 


0.033 


0.040 


0.013 


0.020 


b- likeness of jet 


0.009 


0.009 


0.008 


0.008 


Background normalization 


0.002 


0.004 


0.000 


0.003 


Background shape 


0.035 


0.031 


0.019 


0.013 


Total 


0.057 


0.060 


0.028 


0.029 



a top mass independent value for F + of zero, whereas we 



see a small influence of the top quark mass on our mea- 
surement of F+. For a shift of +5 GeV/c 2 (-5 GeV/c 2 ) 
in the top quark mass we estimate a deviation in F + of 
+0.008 ± 0.003 (-0.008 ± 0.003). 



VIII. RESULTS 

We have presented a method for the measurement of 
the fractions Fo and F + of longitudinally polarized and 
right-handed W bosons in top quark decays using a se- 
lected data sample with an integrated luminosity of ap- 
proximately 318 pb" 1 collected with the CDF II detector. 

Taking the systematic uncertainties into account, as- 
suming a top quark mass of m t — 178 GeV/c 2 , and as- 
suming that the non-measured fraction is equal to the 
standard model expectation, the final result for the frac- 
tions of longitudinally polarized and right-handed W 
bosons is 

F = 0.85±o.22 ( stat ) ± °- 06 ( s y st )> 
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F+ = 0.051^5 (stat) ± 0.03 (syst). 

We obtain an upper limit on the fraction of right- 
handed W bosons of F+ < 0.26 at the 95% CL. The 
systematic uncertainties are incorporated by convoluting 
L(F + ) with a Gaussian with a mean of zero and a width 
equal to the total systematic uncertainty. 

Figure 2a (b) shows the negative log-likelihood as a 
function of Fq where the minimum represents the 

result of the fit. Our method provides the possibility 
to correct the distribution of observed cos 6* for the se- 
lected sample for acceptance and reconstruction effects. 
Figures 2c and 2d show the unfolded distribution, nor- 
malized to the theoretical ti cross section, in comparison 
with theoretical predictions for standard model and a 
V + A model in the case of F + . As one can see, the 
observation is compatible with the standard model pre- 
diction. Also the measured values for F and F + are in 
good agreement with the standard model. 
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